AIP/123-QED 

Method of generating iV-dimensional isochronous nonsingular 

Hamiltonian systems 

A. Durga Devi^, R. Gladwin Pradeep^, V. K. Chandrasekar^ and M. Lakshmanant 
Centre for Nonlinear Dynamics, School of Physics, 
Bharathidasan University, Tiruchirappalli - 620 024, India. 
(Dated: July 20, 2012) 

Abstract 

In this paper we develop a straightforward procedure to construct higher dimensional isochronous 
Hamiltonian systems. We first show that a class of singular Hamiltonian systems obtained through 
the Q-modified procedure is equivalent to constrained Newtonian systems. Even though such 
systems admit isochronous oscillations, they are effectively one degree of freedom systems due to the 
constraints. Then we generalize the procedure in terms of Qj-modified Hamiltonians and identify 
suitable canonically conjugate coordinates such that the constructed Qj-modified Hamiltonian is 
nonsingular and the corresponding Newton's equation of motion is constraint free. The procedure is 
first illustrated for two dimensional systems and subsequently extended to iV-dimensional systems. 
The general solution of these systems are obtained by integrating the underlying equations and is 
shown to admit isochronous as well as amplitude independent quasiperiodic solutions depending 
on the choice of parameters. 

PACS numbers: 
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I. INTRODUCTION 



Nonlinear oscillators model many physical phenomena and the study of such oscillators 
plays an important role in the theory of nonlinear dynamical systems. Various methods 
have been developed to study nonlinear oscillator equations. The dependence of frequency 
on the amplitude of oscillation is generally considered to be a fundamental property of 
a nonlinear oscillator. However, recently several nonlinear oscillators have been identified 
whose frequencies are constants independent of the initial condition, and are known as 
isochronous systems. Moreover, several procedures have been developed to construct and 
identify classes of .sochronous osci.late.rs. In particular Ca.oge r0 a n d Ca.oge r o a n d 

Leyvraz |3H7| have developed many techniques to generate isochronous oscillators and one 
among these is the method of transforming a real autonomous Hamiltonian H(p, q) into 

a el-modified Hamiltonian H = ^(H(p,q) 2 + Q 2 Q(p,q) 2 ), where p = (p±,P2, Pn) and 

q = (gi, 92, ...9at)- Note here that H is the Hamiltonian of the one dimensional harmonic 
oscillator Q + Q 2 Q = 0, and (H, Q) are the canonical conjugate momentum and coordinate, 
respectively. It is interesting to note that the nonlinear equation obtained from this Q 
modified Hamiltonian is isochronous and the system evolves periodically with the fixed 
period T = 2tt/Q. 

Recently Chandrasekar et al proposed a procedure to generate scalar isochronous sys- 
tems recursively from a given Hamiltonian [8j. In a recent paper, a procedure to identify 
several classes of coupled isochronous systems has been proposed and these systems are con- 
structed using certain nonlocal transformations on a system of harmonic oscillators [12| . In 
the present paper, we briefly analyse the above mentioned procedure of el-modified Hamil- 
tonian developed to generate a system of iV-dimensional isochronous oscillators. We find 
that while the systems generated by this procedure may be isochronous, the underlying 
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sequently the Newton's equation of motion becomes constrained with (N — 1) holonomic 



constraints so that the resultant system corresponds effectively to a one degree of freedom 
system. Without these constraints, the Newton's equation of motion will admit unbounded 
solution, while with the introduction of the constraints the solution becomes bounded and 
isochronous. In this paper we show how one can develop a procedure which leads to non- 
singular Hamiltonian systems such that the Hamilton's equation or Newton's equation of 
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motion are constraint free, but which can still admit isochronous oscillations depending on 
the choice of parameters. For other choices one can obtain amplitude independent quasiperi- 
odic solutions. The procedure is first demonstrated for two-degrees of freedom systems and 
subsequently extended to iV-degrees of freedom systems. 

We organize our study as follows. In section 2, we briefly discuss the procedure of Q- 
modified Hamiltonian developed by Calogero and Leyvraz [7|] and note that the deduced 
Hamiltonian is a singular or degenerate one. In section 3, we modify this procedure suitably 
and use it to construct two-degrees of freedom isochronous systems which are free from 
constraints. We obtain the general solution and find it to be isochronous or quasiperiodic 
for appropriate choice of parameters. In section 4, we extend the procedure to iV-degrees of 
freedom to construct isochronous systems. Finally, in section 5 we summarize our results. 
In the Appendix we point out a procedure to obtain the general solution of certain coupled 
nonlinear differential equations which are of interest for the present study. 



II. THE METHOD OF ^-MODIFIED HAMILTONIAN 



Let us consider an iV-dimensional system with a Hamiltonian of the form 



N 



H(p,g) = ^2a n p n [ 



n=l 



dQ(q) 
dq n 



(1) 



where p = (pi,p 2 , ...,Pn) and q = (<Zi, <72> •••> (Ziv) and Q(q) is an arbitrary function of the 
canonical coordinates q' n s, n = 1,2, ...N. The form of the Hamiltonian is chosen in such 



N 

n=l tt n 



1. In other 



a way that the Poisson bracket {H, Q} is unity for the condition 
words, the Hamiltonian H(p,q) and Q(q) are the canonically conjugate variables for the 
choice Y!l=i a n = 1- 

Using the above canonically conjugate coordinate and momentum variables the procedure 
of Calogero and Leyvraz |7J allows one to construct the f2-modified Hamiltonian H(p, q; f2) 
through the relation 



H(p,q-,n) 



1 
2 

1 
2 



{H(p,q)) 2 + Q\Q(q)f 
f fdQ(q)Y l 



Q 2 Q(qf 



(2) 
(3) 
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where f2 is a constant. The above Hamiltonian produces the following 2N first order canon- 
ical equations of motion for the canonical coordinates q n and p n , 

-dQiqy- 1 



Pn 



dq r , 



n 



1,2,. ..,N 



N 



m=l 



i. &rriPm 


IdQ(q)- 


-2 


\d 2 Q{qY 






dq m 




_dq m dq n _ 





dQ(q) 
dq n 



(4a) 
(4b) 



Now, it is known that a given Hamiltonian is known as a singular /degenerate Hamiltonian 
if the determinant of the Hessian matrix is zero, that is the condition 



d 2 H 



dpidpj 



0, i,j 



1,2, 



N, 



(5) 



is satisfied |9l-lllj. Now substituting the form of the ^-modified Hamiltonian H given by 
equation ([3]) into the Hessian matrix, we get 



d 2 H 



dpidpj 



aid?, 



CL\CLN 



dqi J 



dQ 



dQ 
dq N 



a\(l2 



a 2 a 3 



dqi J KdqN J \Qq 2 

One can see that taking out the term Qif^ from each row, where i is the row index, all the 
rows become equal, that is 

-l 
-l 



a 



8Q \ 1 f 9Q N 
9q2 ) \dqi 

2 fdQ^- 2 



2 ^ dq2 

dQ_ \~ ( 9Q N 
dq 2 / \ 9q 3 



dQ 



dQ 
dq n 



CL\(lN 



0>2 a N 



dQ 
dqN 

dQ 
dqN 

a ''' ( '- x I Th, 



a 



-l 



-l 



dQ 
dqi 

dQ 
dqi 

dQ 
dqN 



-1 
-1 
-1 



dQ 



N \ dq~ 



.(6) 



d 2 H 



dpidpj 



N 

ndi 

i=1 \dq t 



dQ 
dqi 

dQ 
dqi 

"Jig J "2 



dQ 

C)(I2 



«1 



dQ 
dqi 



a 2 



dQ 
dq2 

dQ 
<>•!> 



dQ 
dq2 



-1 
-1 



a N 
a N 

a N 



dQ 
dqN 



dQ 
dqN 
dQ 
dqN 



dQ 
dqN 



-1 
-1 
-1 



0. 



(7) 



which implies the Hamiltonian ([3]) is singular /degenerate. Consequently, the Newton's equa- 
tion of motion for such a singular Hamiltonian system is accompanied by a system of con- 
traint equations, that is the corresponding Newtonian equation is a constrained system, as 
we see below. 
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In Ref. [jj, the Newton's equation of motion is obtained by differentiating the first 
equation (Hal) with respect to and replacing H = —Q 2 Q(q). It is given in the form 



dQ(q) 
dq n 



A' 



-q, 



m=l 



dq n q m J 



+ n 2 a n Q(q) = 0,n = 1,2,..., N. 



(8) 



Equivalently, one can obtain (JS} by differentiating (Ha]) with respect to £ and using ( lib]) for 
in it. 

However, one can easily check that not all the coordinates qi, i = 1,2, ... ,N are inde- 
pendent: there are (iV — 1) holonomic constraints existing between them. This can be easily 
seen by the fact that from (Ha|) it follows that 



dQ 



a n \dq : 

Consequently we have the relations 

1 . 9Q{q) 1 . dQ(q) 



q n = H, for each n = 1,2, ... ,N 



—qr 



dqi dj * dqj 

On integration, one obtains a set of (N — 1) functional relations on the coordinates qi, 



0, j=2,3,...,N. 



(9) 



(10) 



i = l,2,...,JV: 



dQ 



dqj^- = Cj, j — 2,3 



,...,N, 



(11) 



where Cj's are constants. Equation (ITT]) obviously constitutes a set of (N — 1) holonomic 
constraints on the coordinates qi. These constraints have to be supplemented with the 
equation of motion (JSJ) which together constitute the Newton's equation of motion. The 
resulting system is then equivalent to the Hamilton's equations @. One may also note 
that equation <^ on its own merit, that is without the constraints (ITT]) , corresponds to the 
nonsingular Hamiltonian 



2 2 
<HPi 



,i=l 



dQ(q) 
dqi 



(12) 



with the associated canonical equations 



qi = ®iPi 

N 



9Q(q) 
dqi 



Pi 



9Q(q) 
dqj 



" 3 «2 



dqjdqi ~ dqi' 



(13a) 
(13b) 
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A. Specific Examples 



In order to illustrate the above points clearly, we consider a specific choice of Q(q) which 
has been considered in Ref. [7|, namely 

N 

Q(q) = Y, b ^ m > ( 14 ) 

m=l 

where 5 m 's are arbitrary real parameters and fc m 's are such that (l/fe m )'s are positive integers. 
Correspondingly the Newton's equation of motion (without the constraints) (jSJ) takes the 
form 

N 

knbn [q k n~ l qn + {K ~ l)^"^] + ^ b m q k ™ = 0. (15) 

m=l 

The associated (N — 1) holonomic constraint conditions ( TTTj) become 
The singular Hamiltonian (J3J) for the choice (|14p becomes 

JV \ 2 / JV \ 2 1 AT 



vm=l / I n=l 



and the canonical equations become 



h ( N \ N 

Err - fi W 4 E^. . = i,2,... ; at. (isb) 

On the other hand, the nonsingular Hamiltonian corresponding to the Newton's equation of 
motion ffl5|) without the constraints ffl6|) is 



H 

2 



with the corresponding Hamilton's equations 



i "i 



9q 2k * ' 



/'' = I^TT " n 2 **^" 1 E i = l,2,...,J\T. (20b) 
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We now demonstrate the isochronous nature of the singular Hamiltonian system (fl8l) or 
equivalently the Newton's equation with constraints f ri5|) -f fT6|) . and the nonisochronous and 
unstable nature of the solutions for the nonsingular Hamiltonian system fll9p or equivalently 
the Newton's equation without the constraints ffTB]) . To be clear we first discuss the N = 2 
degrees of freedom special case and then the N arbitrary case with a simple choice of 
parameters. 



1. N = 2 case 

Let us consider the specific two degrees of freedom Hamiltonian H = + -^ 2 - which is 
canonically conjugate to the collective coordinate Q = (q\ + q%), that is the Poisson bracket 
{H,Q}=1. Here we have chosen the constant parameters a\ = a 2 = 1/2, k\ — k 2 = 1 in 
equation ( TI4")) for simplicity. 

Considering H as the momentum and Q as the canonically conjugate collective coordinate 
and substituting them into the fi-modified Hamiltonian fl3]), we get 

H=\\H 2 + ^Q 2 ] = + + n\ q \ + gin (21) 



The corresponding canonical equations of motion are 

H qiP2 + Q2Pi . H qiP 2 + q2Pi ,,- )0 , 

2i = T~ = Ta 2 » Q2 = — = —— a — , (22a) 

4gi 16gfg 2 4g 2 16gi<?2 

Pi = (P2?i + Pi<? 2 ) - 2ft 2 gi (g 2 + g 2 ) = - 2Q 2 g x Q, (22b) 



(P29i + Pi<? 2 ) - 2n 2 qi (qf + g 2 ) = ^> 2 - 
16gfg 2 4gf 



p 2 = -p-g (pig 2 + P2?i) - 2fi 2 g 2 (g 2 + g 2 2 ) = - 2n 2 g 2 Q. (22c) 

iogig 2 4g 2 



From ( 122]) . we can obtain the system of coupled second order ODEs, 

, + q\ + & ^ 

gi - , qi — • [Z6 

4<?i qi 4g 2 q 2 

The general solution of the above system of second order ODEs is (for details see Appendix 
A) 

<h = ^[2(gi (0) 2 ) + H(0) + g(0)(cos(n*) - 1) + Cm 1 ' 2 , (24a) 

92 = ^[2fe(0) 2 ) + g(0) + Q(0)(cos(fit) - 1) - C(0)f/ 2 , (24b) 

where gi(0), 92(0), H(0), C(0) are the four constants of integration fixed by the initial con- 
ditions and Q(0) = 9i(0) 2 + g2(0) 2 . Obviously the solution fT2"4"|) is unbounded in nature. 



However, from equations fl22ajl . it also follows that 

qm = Q2Q2 (25) 

leading to the constraint 

q!-q 2 2 = I, (26) 

where I is a constant. 

Thus the equivalent Newton's equation of motion corresponding to the Hamiltonian sys- 
tem (122]) is the equation (123]) plus the constraint (T26]) . This means that only one of the two 
variables q± or q 2 is independent. Subject to the constraint (T26]) . the general solution of fl23|) 
can be written in the form 

qi = ^[%i(°) 2 ) + H{0) f m + W)(co8(nt) - (27a) 

q 2 = ^[2fe(0) 2 ) + H(0) f m + Q(0)(cos(nt) - l)] 1 / 2 , (27b) 

where <7i(0) 2 — ^(O) 2 = I, which is periodic, bounded and isochronic, but only with three 
integration constants H(0), q\(0) and 92(0). One can also check that the solution ( l27j) is also 
the solution of the Hamilton's equations ( 1221 . 

Further, we also note that the Newton's equation ( 1231) on its own merit without the 
constraint ( 126]) admits a different Hamiltonian of the form 

" 4 (S + i + ^ + * 2 > 2 ) • < 28 > 

which is not equivalent to the Hamiltonian (|2"T]) . The corresponding canonical equations are 



Pi ■ V2 {oa , 

qi = w q2 = w (29) 



2 2 

Pi = Ps- m qi {qi + g 2 ), p 2 = ^ - 2fig 2 (g! 2 + g 2 2 ), (30) 

8qf 8ql 
which is now equivalent to the Newton's equation ( 123]) alone without the constraint ( 126]) . 
Obviously this unconstrained system admits the unbounded solution ( 1!Mj) . 



2. N-arbitrary case 

In the general case, we consider the singular Hamiltonian ( ITTj) whose canonical equations 
are given by Eqs. ( 1T8]) which lead to (N — 1) constraints given by Eq. ( 1T6]) . The corre- 
sponding Newton's equation of motion is given by Eq. f|T5]) along with the constraints ffTB"]) . 
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Integrating now the Newton's equation (fI5I) without the constraints ({16]) . as pointed out in 
Appendix A, one obtains the general solution as 



q n (t) = q n (0) 1 



1 



b n (g„(0))' 



H(0) 



sin fit + Q(0)(cosfit - 1) + C n (0)t 



(31) 



where q n (0)'s, H(0), and C n (0), ^]^ =1 C„ = 0, n = 1,2,...,N are integration constants 
fixed by the initial condition and Q(0) = J2m=i ^m9m(0) fem - Subject to the (N— 1) constraints 
(J16p . the Newton's equation admits the (N + 1) parameter bounded, isochronous solution 



q n {t) = qn{0)[l + -r 



a,, 



bn (gn(0)) fe " 



^( )™M + Q(o)(oos(nt)-i) 



(32) 



which is also the solution of the Hamilton's equations (fT5|) . Note that the Hamiltonian of 
the constraint-free equation (115j) is 




2 2 



gg(g) 



+ n 2 Q(qf 



(33) 



Our above analysis clearly shows that for the singular Hamiltonian systems ([3]), the equiv- 
alent Newton's equation is a holonomic constrained system (with (N — 1) constraint con- 
ditions) admitting isochronous oscillatory solution as the general solution. Consequently, 
the associated system possesses only one independent coordinate variable. On the other 
hand if one considers the Newton's equation (fT5|) without the constraints, the solution is 
unbounded and the associated Hamiltonian is a nonsingular iV-degrees of freedom system. 
In the next section, we describe a procedure to modify this system such that the new system 
with iV-degrees of freedom admits isochronous oscillations. 



III. SYSTEMATIC METHOD TO CONSTRUCT HIGHER DIMENSIONAL 
ISOCHRONOUS SYSTEMS 

Even though the procedure of fi-modified Hamiltonian is an interesting method to con- 
struct A-dimensional isochronous systems, the fi-modified Hamiltonian is singular (as seen 
in the previous section) and the effective degrees of freedom is only one (that is, the number 
of independent coordinate variables is only one due to the (N — 1) holonomic constraints). 
On inspection for the reason behind this fact, we can easily note that the Q- modified Hamil- 
tonian is defined in such a way that the original Hamiltonian H(p, q), which is a function of 
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FIG. 1: Isochronous oscillation of the system (|23p with the constraint (|26p (a) Time series plot for 
qi (b) Projected (q\ — gi) portrait. 

N canonical momenta p and the N canonical coordinates q, becomes the canonical momen- 
tum of an one dimensional harmonic oscillator. This mapping of a N- dimensional system 
onto an one dimensional harmonic oscillator is the reason for the inability to rewrite the 2N 
first order ODEs as a system of N coupled second order ODEs without constraints. One can 
overcome this problem by suitably redefining the ^-modified Hamiltonian. In this section 
we illustrate this for a class of two-dimensional systems and obtain the general solution for 
a specific example. In addition to the periodic isochronous oscillations, the resultant class 
of two-dimensional systems can also admit amplitude independent quasiperiodic oscillations 
for appropriate choice of parameters. 

Let us define the f2-modified Hamiltonian for a two-dimensional system as 

H = ]^[H l { Pn ,q n ) 2 + H 2 { Pn ,q n f + n\Q 1 {q n ) 2 + n 2 2 Q 2 {q n ) 2 l 71=1,2, (34) 

where H±(p n , q n ), H 2 {jp n , Qn) are rea l Hamiltonians of two arbitrary two-dimensional systems 
and Qi and Q 2 are arbitrary functions canonically conjugate to the Hamiltonians Hi(p n , q n ) 
and H 2 (p n ,q n ), respectively. We have effectively mapped the Hamiltonians of two different 
two-dimensional systems to the Hamiltonian of a two-dimensional harmonic oscillator which 
is now the new fij-modified Hamiltonian. We wish to note here that the new f2j-modified 
Hamiltonian H has two natural frequencies of oscillations, Qi and f^, and depending upon 
the ratio between them we obtain periodic or quasiperiodic oscillations. 

Following the definition of the Hamiltonian H(p, q) in the previous section we define the 
new Hamiltonians Hi{p ni q n ) and H 2 (p n ,q n ) as 
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where Q 



dQ l 



1,2, and a«, 6j are new functions to be determined. From the 



solution of the canonical equations of the two-dimensional harmonic oscillator (134j) we know 
that 



Hi = Acos(ttit + 5 X ), H 2 = Bcos(tt 2 t + 5 2 ), 

A B 
Qx = — sin(fiit + 5i), Q 2 = 77- sin(fi 2 t + $2), 



(36) 



where A,B,8i, 5 2 are arbitrary constants. 

Now, one can easily check that the above Hamiltonians Hi and H 2 are canonically con 
jugate to the collective coordinates Qi(q n ) and Q 2 (q n ), respectively, for the choice 

Q2 qi Qlq 2 



QlqiQ2q 2 

ai = , a 2 



bi 



A ' 

Qlq 2 Q2qi 



A ' 

QlqiQ2q 2 



(37) 
(38) 



A A ' 

where A = Qi qi Q 2q2 — Qi q2 Q2 qi 7^ 0. Substituting the above quantities ( 135|) . (|37|) and ( 138|) 
in the f2-modified Hamiltonian ( 13~4"|) we get 



H 



(p,Q2q 2 - P 2Q2 qi ) 2 , + 



(39) 



A 2 A 2 
Here we note that the above Hamiltonian ( !39l is nonsingular and this can be verified from 
the relation for the Hessian, 



d 2 H d 2 H 



dp 2 dp\dp 2 
d 2 H d 2 H 



A ^ 0. 



(40) 



dp 2 dpi dpi. 

The choice A = Qi qi Q 2q2 — Qi q2 Q2 qi = implies Qi(q) = f(Q2(q)), where / is an arbitrary 
function and for this case the fi-modified Hamiltonian (|34"|) reduces to ([3]), which we have 
already shown to be a singular one. 

The canonical equations of motion corresponding to the above Hamiltonian (139]) are 
1 



qi 



A 2 
1 



(Q2q 2 (PlQ2q 2 ~ frQ^) ~ Q2q 2 (P2Q qi ~ PlQlq 2 )) 



q.2 = ~JT 2 (Ql P2 (P2Ql qi -PlQlq 2 ) - Q2 Ql (PlQ2q 2 ~P2Q2 qi )) 



A 2 
1 



(41a) 
(41b) 



Pi 



A 3 



A(p2<9l 9l ~ PlQlq 2 )(P2Qlq iqi ~ PlQl qi q 2 ) ~ (P2<5l 9l ~ PlQlq 2 ) 
X {QlqiQ2 qi q 2 + Q2q 2 Qlq 1 q 1 ~ Qlq 2 Q2q iqi ~ Q2q 1 Qlq 1 q 2 } ~ (PlQ2q 2 ~ P2Q2 qi ) 

xA(piQ 2qiq2 — p 2 Q qiqi ) — ij>iQ 2q2 — P2Q2 qi ) 2 {Ql qi Q2 qiq2 + Q2 q2 Ql qiqi 



Qi q2 Q2 q i q i Q,2 q \Q\ q \ q2 } 



+ QiQiQi qi + Q 2 Q2Q2 qi -, 



(41c) 
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and 



P2 



1 

A3 



A(p 2 Ql qi ~ PlQlq 2 )(p2Qlq iq2 -PlQlq 2 q 2 ) ~ (P2<2l 9l ~ PlQlq 2 ] 



X (QlqiQ2q 2 q2 + Q2q 2 Ql qi q 2 ~~ Q\q 2 Q2q\q 2 ~ Q2q 1 Qlq 2 q2) ~ (PlQ2q 2 ~P2Q2qx, 

xA(p 1 Q 2q2q2 — P2Q2 qi q 2 ) ~ (PlQ2q 2 ~ P2Q2 qi Y {Ql qi Q2q 2 q 2 + Q2q 2 Qlq 1(l 

l2. 



192 



~Qlq 2 Q2qiq 2 Qlq 2 q 2 Q2q\\ 



tiiQiQiq 2 + n 2 2 Q 2 Q 2q2 . (4id) 

Rewriting the above Hamilton's equations of motion we get the following system of two 
coupled second order ODEs, 



Qi 



\fc=l j=l J 



1,2, 



(42) 



where 



-1 



Qlq i+1 Q2q i+1 
Qlqjqk Q2qjq k 



B; 



:-ir 



Qj = Qj(qi,q 2 ),2 + i 



(43) 



Qi&i Q 2 Q 2 

Qlq i+ l Q2q i+1 

In order to obtain the explicit general solution of ( 141]) or fj42l) one has to fix the form of 
Qi and Q 2 in the above equation such that the resultant solutions are analytic and single 
valued. For illustration we choose 



Qi = hq'l 1 + k 2 q 2 2 Q 2 = hq r 1 1 +hq?, 



(44) 



where r% and r 2 are such that (1/ri) and {l/r 2 ) are positive integers (see below) so that the 
resultant solution is single valued and analytic. One can also choose other forms of Qi(q n ) 
and Q 2 (q n ) as well but we stick to this form for the sake of simplicity. Substituting the 
above forms in fl39l) we get 



6=1 

2 



{ktf>\r 2 q x ri - hp 2 nq 2 



l-r 2 \ 



(k 2 p x r 2 q\ Tl - k x p 2 r x q 2 



1— r 2 ^ 



92r x r 2 

•\2fi„ „n i h. „r 2 \2 , C)2(i, n i r, „r 2 \2 



92T\T 2 



+ni(k 1 q r 1 1 + k 2 q r 2 2 Y + nHhq? + k 4 q? 
where g 2 = k 2 k 3 — k\k± ^ 0. The corresponding equations of motion are given as 
k 2 k 4 (Ql - Q 2 2 )ql- ri q r 2 2 



(45) 



Qi 



ri92 

+qi X {{hhttl - k 2 k 3 n 2 2 )ql + (-k 2 k 3 + kxh)(-l + r x )r iq l} 
1 



(46a) 



qi 



r292 



hh^i-nDqi^ 



+q 2 1 {(-k 2 k 3 Vt\ + kikiVt 2 2 )ql - (k 2 k 3 - kik A )(-l + r 2 )r 2 q 2 ] 



(46b) 



12 



MOO 



40 



son 




-O.l 1 — 
0.4 



31 



(c) 



FIG. 2: Isochronous oscillations admitted by Eq. (|46p for the choice ill : = 1 : 1, r i= r 2 = \ 
(a) Time series plot for q\ (b) Projected (gi — <ji) portrait (c) Poincare section 



The general solution of the above system of equations can be obtained by either explicitly 
integrating the Hamilton's equations of motion or by substituting the relation ( )36|) in ( )44|) 
and solving the resultant algebraic equations for q\ and q 2 . The explicit general solution of 
[6]) can be given as 



Qi 



Q-2 



sin(^ 2 t + 8 2 ) - sin^t + 8 1 



jfc sin(^ + sin(n 2 t + 5 2 ) 



, {k 2 k 3 - kik 4 ) ^ 0, 



(47a) 



(47b) 



k 2 k 3 - kik± 

We note here that the above solution contains the required number (four) of arbitrary 
constants, namely A, B, 5i, and S 2 . The obtained solution ( 147|) is analytic and bounded and 
exhibits oscillatory behaviour for the choice 1/ri and l/r 2 are positive integers. 

In Figs. 2 and 3, we plot the solution of Eq. (T4"rJj) for two different ratios of Qi and 
Q 2 with the parametric choice r\— r 2 = |, ki=k 2 =k^=l and k 3 = 2. In Fig. 2, we choose 
Qi : Q 2 = 1 : 1 which results in isochronous periodic oscillations. In Fig. 3 we choose 
Hi : £l 2 = I : V% which results in quasiperiodic oscillations with amplitude independent 
frequencies. 
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400 



(a) 




FIG. 3: Quasiperiodic oscillations admitted by Eq. ([4"6|) for the choice f2i : O2 = 1 : v2, T\= 
T2 = 5 ( a ) Time series plot for gi (b) Projected (gi — q\) portrait (c) Poincare section 

Having developed the method of constructing isochronous two dimensional oscillators, we 
extend the procedure to iV-dimensions in the next section. 



IV. N-DEGREES OF FREEDOM GENERALIZATION OF HAMILTONIANS 

One can generalize the procedure of constructing isochronous Hamiltonian systems to 
N degrees of freedom. In order to do so, we have first extended the procedure to three 
dimensions and then generalized it to N dimensions. For the sake of brevity we do not 
present the results of three dimensional systems here. Proceeding further, to construct the 
N-dimensional isochronous Hamiltonians let us consider the following Hamiltonians 

H i = J^ a ijPj[-^J , i = l,2,...,iV (48) 

which are canonically conjugate to the corresponding collective coordinates Qi, i = 
1,2, ... N. Here the functions a^-'s are given by 
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where 



Qij Q 



Q(i+l)q j+1 Q(l+ 
Q(i+2)q j+1 Q(i+2)q J+2 



QNq j+1 QNq 



J+2 



i+l)q N i+l) qi Q(i+l)q 2 ■ ■ ■ Q(i+l)qj-i 
i+2)q N Q(i+2) qi Q(i+2)q 2 ■ ■ ■ Q(i+2)q j - 1 



QNq N QNqi QNq2 



QNq,--, 



(49) 



Q 



dQi 

dq 3 



A 



Qlqi Qlq2 ■ ■ ■ QlqN 
Q2 qi Q2q 2 ■ ■ ■ Q2q N 

QNq, QNq 2 ■ ■ ■ QNq N 



1,2,. .N, 



(50) 



with N + k 



k — ; 1. 2. 



N — 1 . Here Qij isaniV — lxiV — 1 determinant and A is 
an N x N determinant. In order to construct the iV-dimensional isochronous Hamiltonian 
system, let us substitute the H^s given by the relation (148]) into the Hamiltonian of the 
iV-dimensional harmonic oscillator specified by 

N 

(51) 



1 N 



i=l 



where Hi and Qi 



1,2,..N, are the canonical momenta and coordinates, respectively 



Using this Hamiltonian H and the corresponding Hamilton's equations of motion we find 
the following system of iV coupled second order ordinary differential equations (ODEs), 



qi 



l / N N \ 

^ ( Yl A v>°Wk + B { ) , 

U=l ]=! / 



1,2, 



, N, N > 2 



(52) 
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where A^ and By, j, k = 1, 2, . . . , N are determinants of the form 



.4 



ijk 



B, = - 



Qlq i+1 Q2q i+1 

Qlq i+ 2 Q2q i+2 

QlqN Q^qN 

Qlqi Q,2q x 



Qlq t -! Q2q i ^ 1 
Qlqjq k Qlqjqk 

Q\Vt\ Q2&2 
Qlq i+1 Q2q i+1 
Qlq i+2 Q2q i+2 

Qlq N Q2q 

N 

Qlqi-1 Q2q i _ 1 



■ ■ QNq i+1 

■ ■ QNq i+2 

■ ■ QNq N 

■ ■ QNqi 



■ ■ QNq,.! 

■ ■ QN qj q k 

Qn Qi+1 

QNq i+2 



■ ■ QNq N 

■ ■ QNq,.! 



(53) 



NxN 



(54) 



NxN 



The general solution of ( 1521) can be found after choosing the forms of Qi(qi, (72, • • • , <Zat) 
and following the procedure discussed for two coupled second order ODEs f H6|) . Using the 
relation 



Ha = A- cos 



and 



Qi 



Sill 



Qit + Si 



where Ai and Si, i = 1, 2, . . . , N, are arbitrary constants, or integrating the Hamilton's 
equations of motion corresponding to the Hamiltonian (l5Tj) . one can easily see that the 
canonical variables qi, i = 1,2, ...,N, also evolve periodically with a fixed period T when 
Qj's are commensurate, for appropriate forms of Qi(qi,q2, ■ ■ ■ ,qN) such that the resultant 
solutions qiS are analytic and single valued. 
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V. CONCLUSION 



In this paper, we have shown that a class of singular Hamiltonian systems obtained 
through the f2-modified procedure are equivalent to constrained Newtonian systems. We 
have made use of the idea behind the procedure of f2-modified Hamiltonian and introduced 
suitable modifications to the procedure and developed a systematic procedure to construct 
a class of nonsingular N- dimensional isochronous systems exhibiting bounded isochronous 
oscillatory motion. The Hamiltonian constructed by this procedure is nonsingular and the 
generated system is free from constraints and the Hamilton's equations and Newton's equa- 
tion of motion are equivalent. The procedure is first developed for the case of two-degrees 
of freedom systems and subsequently generalized to A-degrees of freedom systems. The 
systems obtained through this procedure can also exhibit quasiperiodic oscillations with 
amplitude independent frequencies when the original linear harmonic oscillator frequencies 
are incommensurate. 
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VI. APPENDIX 

The general solution of the system of two-coupled second order nonlinear ODEs (1231) can 
be found by the following procedure. Using the point transformation 



(55) 




(56) 



From the above system of ODEs, one can rewrite 



x + y 



n 2 



(x + y), x-y = 0. 



(57) 
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Defining z\ = x + y and z 2 = x — y we get 

z\ = -fi 2 zi, z 2 = 0, 



(58) 



whose general solutions are 

zi = A sin Qt + B cos fit, z 2 = I\t + I 2 , (59) 

where A, 5, Ji, J 2 are integration constants. Solving for x and y we get 

x = - (A sin fit + 5 cos fit + Jit + J 2 ) , (60) 

y = - (A sin fit + 5 cos fit - lit - J 2 ) . (61) 

Redefining the above arbitrary integration constants as A = ^rr~, B = qiifi) 2 + <? 2 (0) 2 , 
h = C(0), J 2 = gi(0) 2 - g 2 (0) 2 , we get the general solution of Eq. flUD 

i 

?i = ^ (2 gi (0) 2 + H ^ Qt + Q(0)(cosfit - 1) + C(0)t\ 2 , (62) 

g 2 = (2g 2 (0) 2 + + g(0)(cosfit - l) - C{0)t) 2 , (63) 

where gi(0), g 2 (0), H(Q) and C(0) are arbitrary constants and Q(0) = gi(0) 2 + g 2 (0) 2 . The 
procedure can be straightforwardly extended to deduce the general solution of the iV-coupled 
ODEs ©. 

Using the transformation Xj = 6^' we find Eq. ([8]) reduces to the form 

fi 2 * 

m=l 

From the above system of ODEs, one can rewrite 

N N 

} u x m = -fi 2 ^2 Xm > xx-Xj = 0, j = 2, 3, . . . , N (65) 

m=l m=l 

Defining the variables Z\ = J2m=i x m an d z 2j = x\ — Xj, j = 2, 3, . . . , N, we get 

'i\ = -fi 2 ^i, z 2 i = 0. (66) 

Solving the above system of linear ODEs and rewriting the resultant solution in terms of 
the original variables g^'s we get the general solution as 



q n {t) = q n (0) ( 1 + yj^yr f sin fit + Q(0)(cosfit - 1) ) + C n (0)t ) " , n = 1, 2, . . . , if67) 
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where q n (0)'s, H(0) and C n (0), n — 1, 2, . . . , TV are integration constants fixed by the initial 
condition, Eli C n = and Q(0) = E™=i Mm(0) fcm . 
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